We present a molecular dynamics simulation of shock waves propagating in dense deuterium with the electron force field method [JT Su and WA Goddard, Phys. Rev. Lett. 99, 185003 (2007)], which explicitly takes the excitation of electrons into consideration. Nonequilibrium features associated with the excitation of electrons are systematically investigated. We show that chemical bonds in D 2 molecules lead to a more complicated shock wave structure near the shock front, compared with the results of classical molecular dynamics simulation. Charge separation can bring about accumulation of net charges on large scales, instead of the formation of a localized dipole layer, which might cause extra energy for the shock wave to propagate. In addition, the simulations also display that molecular dissociation at the shock front is the major factor that accounts for the “bump” structure in the principal Hugoniot. These results could help to build a more realistic picture of shock wave propagation in fuel materials commonly used in the inertial confinement fusion.